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LOW DIMENSIONAL BORN-INFELD EQUATIONS COUPLED WITH A 
COLLISIONLESS MATTER MODEL 

HO LEE 


Abstract. We consider the Born-Infeld nonlinear electromagnetic field equations and study its 
Cauchy problem in the case that the Vlasov equation is considered as a matter model. In the 
present paper, the Vlasov equation is considered on the so-called one and one-half dimensional 
phase space, and in consequence the Born-Infeld equations are reduced to a quasilinear hyperbolic 
system with two unknowns. A transformation is introduced in order to make the field equations 
easy to handle, and suitable assumptions are made on initial data so that the nonlinearity of the 
field is controlled. 


1. Introduction 

The Born-Infeld (BI) electromagnetic theory [2] was originally proposed as a nonlinear correction 
of the Maxwell theory in order to overcome the problem of infinities in the classical electrodynamics 
of point particles. The underlying idea was to simply modify the classical theory not to have physical 
quantities of infinities, that is the principle of finiteness. It was to replace the original Lagrangian 
density for the Maxwell electrodynamics with a square root form with a parameter 5, by which the 
finiteness of electric fields is ensured. This was the same with the way the special relativity has 
taken into account the finiteness of the speed of light c, i.e., —mc^-y/l — jt? replaced \mv^. The 
exact form of the BI Lagrangian density appears in the next section. 

This theory in recent years again has received much attention since the string theory found its 
relevance to the BI theory [211 US] : the dynamics of electromagnetic fields on D-branes is described 
by the BI theory, and finiteness of electric fields is naturally observed. Be that as it may, this paper 
is not going to the string theory, but the BI theory will be considered rather as a nonlinear version 
of the Maxwell theory as it was the original standpoint of Born and Infeld. Hence, it makes sense to 
choose the Vlasov equation as a matter model since this equation describes the dynamics of classical 
particles. 

In the present paper, we study the BI equations coupled with the Vlasov equation as a matter 
model, and this induces inhomogeneities to the BI equations. For homogeneous cases, we can find 
several results on the Cauchy problem. In [5] , Chae and Huh proved the existence of global classical 
solutions for small initial data. They considered the BI equations as quasilinear wave equations 
and crucially used the null form structure of the nonlinear terms. On the other hand, Brenier |3] 
considered the BI equations as a system of hyperbolic conservation laws. He supplemented additional 
conservation laws to obtain an augmented BI system and discussed conditions for existence of global 
solutions in one dimensional case. The similar frameworks can be found in 12111521123123111 ], and 
for different approaches we refer to [TTl HH US] ■ In this paper, the BI equations will be considered 
as a system of hyperbolic equations as in [3] , but it differs in that inhomogeneous terms are taken 
into account. 

On the other hand, the Vlasov equation has been widely used to describe matters in connection 
with other field equations: electrostatic or electromagnetic fields, non-relativistic or fully relativistic 
gravitational fields, and so on. For detailed discussions on the Vlasov-type equations, we refer to 
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[HIE] and the references therein. The analysis of the Vlasov-Maxwell system is applied to our case, 
and as a simple but nontrivial case we consider the “one and one-half dimensional” case as in M- 
The Vlasov equation in this dimensional case is introduced in Section 2. 

Before we proceed, we briefly discuss some points of this paper. 

1. To make a coupled system of the BI and the Vlasov equations, we should determine an 
analogous equation corresponding to the Lorentz equation in the Maxwell case. In other words, 
we have to derive an equation of motions of particles interacting with the BI electromagnetic field. 
However, as it was pointed out in [5], it is impossible to derive separate equations for the charged 
particles due to the nonlinearity of the field. Instead, we use the “generalized Lorentz force”, which 
is available for small and almost constant helds im, and this will be discussed in Section 2. 

2. As we mentioned above, the BI equations will be considered as a system of hyperbolic equations. 
To better understand its hyperbolic structure, we introduce a transformation in Section 3. This 
transformation is quite tricky but simple. By this transformation, the linear degeneracy of the BI 
system is then easily seen, and moreover it makes the system easy to handle, for instance the proof 
of Lemma 5.6 would be much more complicated without the transformation. 

3. The main difference from the Maxwell case is that we have to control the characteristics of the 
field equations as well as of the Vlasov equation. In Glassey-Strauss’ result [E], spatial and temporal 
derivatives were split into linear combinations of the Maxwell and the Vlasov characteristics called 
T and S derivatives respectively. This was possible because the Maxwell field always propagates 
with the speed of light, while the Vlasov particles do not move with the speed of light. However, it 
will be shown in Section 3 that the BI field propagates with speed — cos/3 or cos a, which may have 
values less than the speed of light c = 1. Hence, we should make a suitable assumption on initial 
data in order to apply the argument of [E], and this will be discussed in Section 4. 

4. In contrast to Glassey-Schaeffer’s result El: we obtain a local-in-time result. In that paper, 
the Maxwell electromagnetic field is estimated as |i3|-|-|i3| <c<oo, (see Lemma 1 and Gorollary 
1 in [14]). and this estimate implies that any particles considered have momentum v with growth 
order v ^ t. However, this is not enough to be applied to our case. If a momentum grows, and 
consequently tends to infinity, then its corresponding velocity will tend to the speed of light, and 
finally a resonance between the particle velocity and the BI field propagation may appear in a finite 
time. This will be discussed in Section 5.2. It seems that more qualitative analysis for the BI and 
the Vlasov characteristics is required to obtain a global-in-time result. 

Notations. We collect some notations which are used in this paper. 

• The speed of light c, and mass m and charge q of each particle are assumed to be unity: 

c = m = q = 1. 

• Greek indices run from 0 to 3 and Latin indices from 1 to 3. The indices are raised or 
lowered by multiplication with the Minkowski metric rjoip = diag(—1,1,1,1), and we use the 
Einstein summation convention such as 

Xay“ = x“ya = ??a,sx“y^ = -x°y° -h x^^ + x^y^ -b x^y^. 

This notation will appear only in Section 2.1. 

• For a scalar or vector valued function v S d = 1, 2, 3, we define 

V 

y/TT^- 

Then for a momentum G K^, its corresponding velocity is defined by v. 

• The usual norms, 1 < p < oo, are used. For a function f = f (x) defined on 

X = (xi, • • • , Xr) G X • • • X 
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its norms are defined as follows: for 1 < p < oo and each i = 1, - ■ ■ , r, 

l|f||p = / |f(x)|P(ixi • • • dXr, 

X-'-XW^r 

||f(xi)||P= / \f{x)\P dxi ■ ■ ■ dxi^idx^+i ■ ■ ■ dxr, 

and similarly for p = oo and each i = 1, - ■ ■ , r, 

||f||oo =sup{|f(y)| : y e x ••• x M"*’'}, 
l|f(xi)lloo =sup{|f(y)| lyeM'^i x ••• xK‘^'',y, = xj, 

and so on. 

This paper is organized as follows: we first set up the Cauchy problem for the coupled system in 
Section 2, where the Born-Infeld equations will be properly coupled to the Vlasov equation in the 
one and one-half dimensional case. In Section 3, we introduce a transformation which transforms 
the Born-Infeld equations into a quasilinear system in a diagonal form. In Section 4, we state the 
main result of the present paper. Suitable assumptions will be made on initial data in this section. 
Section 5 and 6 are devoted to the proof of the main theorem. 

2. Problem setting: coupled system in one and one-half dimensions 

In this section, we study the Born-Infeld equations and the Vlasov equation to consider their 
coupled system in a low dimensional case which is called the one and one-half dimensions. In the 
first part we briefly review the Born-Infeld electromagnetic theory, and then in the second part we 
consider the one and one-half dimensional case for the coupled system. 

2.1. The Born-Infeld equations. One can obtain field equations by constructing a suitable La- 
grangian density C{F) and then applying the Euler-Lagrange equations to it. In the case that 
charged particles are given as source, the interaction between the particles and field must be consid¬ 
ered in the construction of a Lagrangian density. In this paper, we simply construct a Lagrangian 
density by adding the standard interaction term Aaj°‘ to the Lagrangian as follows: 

^total 

where C{F) is the Lagrangian density for the field itself, and Aa is a four-potential from which the 
electromagnetic field tensor is defined by E = Fap = daA/^ — d/^Aa- In the Born-Infeld theory, the 
Lagrangian density jC(F) is defined as follows: 

£(F) = 6^ det ^Pa/3 + 'j , 

where r/a /3 = diag(—1,1,1,1) is the Minkowski metric, and 6 > 0 is a parameter which measures the 
nonlinearity of the field. For simplicity, we take & = 1 in the present paper. The Euler-Lagrange 
equations now give the equations of field: 

= r, where 

To have explicit formulae for the field, instead of the tensor form, we use T), and H in place 

of Fajs and Ga/s as follows: 

3 3 

^ij ^ ^ ^ijk^k-) G^q and Gij ^ ^ 
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where e^-fc is the completely antisymmetric tensor such that £123 = 1, £213 = —1, and so on. Then, 
we get the following system of PDEs: 

V X E = -dtB, V • B = 0, 

12 11 t , , 

^ ’ V-D = p, VxH = dtD+j, 


where p = and j = are the charge density and the current density respectively. Note 

that we obtain Maxwell’s equations by setting D = E and H = B. The main difference from the 
Maxwell theory is the following nonlinear constitutive relations between E, B, D, and iJ: 

^ _dC E+{E-B)B 

~ BE ~ - \E\^ + |B|2 - {E ■ BY' 

^ _ dC _ B-{E- B)E 

dB~ \/l- \E\^ + \B\^ - {E ■ 

In the above relations, D and E[ are given as functions of E and B. However, we convert them into 
a more convenient form for later use, that is, we rewrite E and H as functions of D and B as follows: 




D-{Dx B) X B 
+ \D\^ + |H|2 + \Dx 
B + {D X B) X D 
^/l + \D\^ + \B\^ + \Dx B\^’ 

which can be verified by direct calculations (see Chapter 20 of [H]). As a result, we obtain the field 
equations (2.1|-(2.3) for given charge and current density. 


( 2 . 2 ) 

E = - 

\ 

(2.3) 

H - 


We have now obtained the equations of field for given sources p and j. On the other hand, in order 
to make a self-consistent coupled system, we have to find the equations of motions of particles for 
given fields D and B. In the linear electromagnetic theory, it is well-known that particle trajectories 
are determined by the Lorentz equation. Let p be a momentum of a particle with charge q and v 
the corresponding velocity. Then, the Lorentz force exerted on the particle is given as follows: 

dp p 

-rr = q{E + v X B), where v = - 

dt -x/i + IpP 

However, it is not easy to derive explicit equations of motions of particles when they are interacting 
with the Born-Infeld electromagnetic field due to the nonlinearity. We refer to (3 03 El for alternative 
methods. In the case that a small electromagnetic field is considered, we can use the “generalized 
Lorentz force” as in 017], where the author studied the dynamics of a particle “dyon” interacting 
with the Born-Infeld field. Dyon is a hypothetical particle which has both electric and magnetic 
charges. Let Qe be electric charge and qm magnetic charge of a particle. Then we have the following 
generalized Lorentz force: 


^ = qe{D + v X B) + qm{B -vxD). 

In the present paper, we will not consider the particles having magnetic charges, hence we take the 
equation of motions of particles as the following equation: 


dp 

dt 


(2.4) 


qe{D + V X B). 
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2.2. One and one-half dimensional case. We now set up the problem for the one and one-half 
dimensional case by following the framework of [T3], where the authors studied the Vlasov-Maxwell 
system. In this case, the spatial variable is a; = (x, 0, 0) G and the momentum variables are 
V = (vi,V2,0) G Hence, the electric field U and the magnetic field B are given by 

H = (Di, iD 2 ,0) € and B = (0,0,B) G 

where Bi, H : / x i—)■ i = 1,2, for a suitable time interval / C [0, oo), and the charge density 

p = p{t, x) and the current density j = j{t, x) are given by 


pe 


and j = (ji,j 2 , 0 )e 


Then, the field equations (2.11 are reduced to 

dtDi = -ji, 

dtD2 = -d,^H - j2, 

dxDi = p, 
dtB = -dxE2, 

and we note that Di is obtained from the first and the third equations: 

Di{t,x) = - ji{s,x)ds or Di{t,x) = / p{t,y)dy. 

Jo J — oo 

These are consistent with each other because the densities p and j will be induced by the Vlasov 
equation, from which the continuity equation “dtp + dxji = 0” is naturally obtained, which makes 
the consistency between them. We will use the latter formula to define Di, and now we have only 


two equations of helds. The constitutive relations (2.2|-(2.3) are reduced to 

{l + B^)D 


E = 


H = 


^1 + \D\^+B^ + \D\^B^ 

o + ms 

^1 -f \D\2 + B^ + \D\^B^ 




D 






B 


Vi + be 


We now consider the Vlasov equation. In contrast with El: where the Lorentz force was used, we 
take the generalized Lorentz force (|2.4[) to have the following equation: 


dtf + vidxf + (Di + V2B, D2 - viB) ■ V„/ = 0, 

where v denotes the momentum of a particle, and v is the corresponding velocity: v = u/(l-|-|wp)^/^. 
The charge density p and the current density j are induced by / as follows: 

p{t,x)= / f{t,x,v)dv — n{x) and j{t,x) = / vf{t,x,v)dv, 

Jr2 

where n{x) is a neutralizing background density in the sense that 

[ [ f™{x,v)dvdx= f n{x)dx, 

Jr Jr^ Jr 

where /*" is an initial data of /, and the above quantity is preserved in time by the mass conservation 
of the Vlasov equation. 

As a result, we will study the Cauchy problem of the following system of PDEs. For f > 0, x € K^, 
and V G K^, consider the distribution function / = f{t,x,v), the electric field D — D{t,x) G K^, 
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and the magnetic field B = B{t,x) G satisfying 

(2.5) dtf + vM + (Bi + V 2 B, D 2 - viB) ■ V,/ = 0, 


( 2 . 6 ) 


(2.7) 


dtD2 + dx 


dfB + dx 


V^l + \D\^ + B^ + \D\^B^ 
{1+B'^)D2 


= -J2, 


= 0 , 


yi + IDP +^2 + |£)|2S2 

where the charge density p = p(t, x) and the current density j = j{t, x) are given by 

(2.8) p{t,x)= / f{t,x,v)dv — n{x) and j{t,x)= / vf{t,x,v)dv, 
and the first component of D field is defined by p as follows: 

(2.9) Di{t,x)=l' p{t,y)dy. 

J —00 


Remark 2.1. The system (2.5)-(2.9) can be thought of as a nonlinear version of the one and one- 
half dimensional Vlasov-Maxwell system m. Maxwell’s equations have been replaced by a system 
of hyperbolic equations (2.6)-(2.7). 


3. Transformed system 


In this section, we study the field equations (2.6)-(2.7) for given charge density p and current 


density j. We introduce a useful transformation by which we can analyze (2.6)-(2.7) more effectively. 


Note that (2.6)-(2.7) is a hyperbolic system with two unknowns D 2 and B for given p and j since 


p induces the first component of Zl by (2.9). 


3.1. Transformation of the field equations. As a first step, we rewrite the hyperbolic system 


(2.6)-(2.7) into a quasilinear form. 


Lemma 3.1. Let D 2 = D 2 {t,x) and B — B{t,x) be -solutions to the hyperbolic system (2.6)-(2.7) 
for given -functions p and j with Di defined by (2.9). Then, the hyperbolic system is rewritten as 
the following quasilinear form: 


(3.1) 


dt 


D 2 

B 


An Ai 


dx 


A 21 A 22 

where the components of A= {Aij)ij^i ^2 are given by 

D 2 B 


D 2 

B 


Cip- j2 

C 2 P 


Aii = 


yi + |L»|Vl + B2’ 


Ai 2 = v'l + lzip 


A 21 = Vl + R2 


(1 


Dl) 


( 1 + |LI|2)3/2’ 


A 22 — 


1 

(1 + R2)3/2’ 

DoB 


and Ci, i = 1,2, are given by 
Ci=- 


DiB 


Vl + |Z?|Vl + R2 


and C 2 = \/l + i?2 


Vl + |Z?|Vl + R2 ’ 


D^ D 2 


( 1 + |D|2)3/2- 


Proof. The proof is an elementary calculation, so we only remark the fact that ^-derivatives of Di 
have been replaced by p. Hence, we obtain a quasilinear form with respect to D 2 and B for given p 
and j. □ 


We next compute the eigenvectors and the eigenvalues of the matrix A. 
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Lemma 3.2. Consider the matrix A in (3.11; 

D 2 B 

^1 + izipvi'+ 52 


/ 


A = 




(l + B2)3/2 


\ 




D 2 B 


(1+11)12)3/2 

Then, its left eigenvectors can be chosen by 

( yr+m -1 \ 

^ 1+|Z1|2 ’ 1+^2^ 

and the corresponding eigenvalues are given by 

D 2 B -y/l + Dj 


and I 2 = 


\/r+^ 1 


Ai — 


^ 1 + |i^| 2 Vl + B 2 


and X 2 = 


l + |Li |2 ’1 + 52 


D 2 B + v^l + £>2 

v/i + |i:>|Vi + B2’ 


Proof. This lemma is proved by direct calculations, so we skip the proof. 


□ 


Remark 3.1. If the fields D and B do not blow up at finite time, then the eigenvectors li and I 2 
are linearly independent, and X 2 is always strictly greater than Ai. Thus, we can see that the system 
([33 is strictly hyperbolic as long as its solution exists. 


Note that (3.11 is a quasilinear hyperbolic system with two unknown functions D 2 and B, and 
this system is transformed into a diagonal form by the following transformation $1 and ^ 2 ' 


$1 : 


/ TT 7r\ / TT 7r\ 

b2-2Jn-2'2) 


$9 


/ TT 7r\ 

V 2’ 2/ 


such that they are defined as follows: for a; = (xi, 0 : 2 ) € and y G we define 


$ii(xi,X 2 ) := arcsin 


+ xf + X 2 


and $ 2 ( 1 /) := arcsin 


yi+ 


yZ 


for 7 = 1,2, where $1 = ($ 11 , $ 12 ). Note that <i)i and $2 are smooth mappings and have smooth in¬ 
verses on their images, $i(IR2) = {((/)i,f) 2 ) : |^i|-|-|^ 2 | < 77 / 2 } and$ 2 (R) = (~77/2,7r/2) respectively. 
We now define the transformed variables as follows: 

6*1 := $ii(Ili, II 2 )) ^2 := <S^i 2 (T>i, D 2 ), and 0 b'-=^ 2 {B). 

In other word, the following relations hold: 

Di . „ D 2 . „ B 


(3.2) 


sin 01 = — , sin 02 = — , sin 0 H = ,_, 

VT+W + Vl + ^ 


which are illustrated in Figure 

On the other hand, since Di is a given function of p, we get the inverses as follows: 


D 2 = -^1 + Df tan 02 and B = tan0B. 


Using the above transformation, we will rewrite (3.1) as a hyperbolic system with two unknown 
functions 02 and 9b- To do that, 0i is regarded as a function of 62 and Di, i.e., the following 
identities will be used in the proof of Lemma 3.3: 


(3.3) 


Di 11 

sin 01 = — , cos 02 and — , =— , cos 02 . 

./i+W 


By the following lemma, the quasilinear hyperbolic system (3.1) is reduced to a diagonal form. 
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Figure 1. The unknowns —7r/2 < 9i < 7r/2, * S {1,2, B}, are well-defined for finite 
Z?i, D 2 , and B. Note that if |0i| -|- |02| —>■ 7i/2, then either Di, D 2 , or both of them 
tend to infinity, which is also the case for 9b- Hence, we can see that 9i ^ ±7r/2 
unless the fields blow up. 


Lemma 3.3. Let D 2 and B be -solutions to the quasilinear hyperbolic system (3.11, and suppose 
that p and j are given -functions satisfying 


(3.4) 


dtDi = -ji and d^Di = p 


for Di defined by (2.9). Consider the transformation (3.2). Then, the quasilinear system (3.1) 
is transformed to the following inhomogeneous linearly degenerate system of a diagonal form for 
a = a(t, x) and (3 = j5{t, x): 


(3.5) dt 


— COS/ 
0 


0 

cos a 


dx 


cos02(fco sin^B 
cos02(fco sin^B 


ki sin 02 + ^2 cos 02 ) 
ki sin 02 -I- k 2 cos 02) 


where a = 92 — 9b and /3 = 02 -I- 9b, and ki, i = 0,1, 2, are given by 

pDi jiDi j2 


kn = - 


i + nr 


ki = 


1 + Df' 


k 2 = - 




Proof. We first transform the eigenvalues A^, i = 1, 2, as follows: 

Ai = sin02 sin0B — cos 02 cos 0b = — cos(02 -I- 0b) = — cos/3, 
A 2 = sin 02 sin0B -I- cos 02 cos 0b = cos(02 — 0b) = cos a. 


We now consider the eigenvectors h, i = 1,2. From the well-known transformation of polar co¬ 
ordinates, {x,y) G <->■ (r, 0) G [0, 00 ) x [0,27r), where x = rcos0 and y = rsin0, we know 
that 

dx9 = — sin0 and dy9 = - cos0. 
r r 

Using the above relations, we first calculate the following quantities. As we can see from Figure 
02 can be thought of a function of Df and II 2 . 

902 ^ 902 9(v/rT)^) 

dDi d{y/l -b Dl) dDi 
1 ■ a 

=- , - Sin 6/9 — , =- 

^/i + W \/i+^ 1 


(3.6) 


cos 02 sin 02 , 
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where we used (3.3). 

(3.7) 

(3.8) 


862 _ 1 

863 1 

Im 


- cos O 2 = 


cos9b = 


\/ r +^ 

1 


VI+ 52 

Hence, the eigenvectors k, i = 1, 2, are transformed to 

f 802 89 b 


1 + H 2 ' 


h = 


\8D2' 8B 


\ , , [ 802 d9B\ 


We now multiply the transformed eigenvector li to (3.1) and obtain 

/ 802 8D2 80b 8B^ 


\8D2 8t 8B 8t 

802 DiB 


+ Ai 


8 D 2 + + 8 D 2 

and then we add the following quantity on both sides: 


802 8D2 _ 89b 8 B 
8 D 2 8 x dB 8 x 

^^2 . d 0 B ^ 
P- ^7=r.?2 - ttttV 1 + 52 


D 1 D 2 


8B 


(1 + |D|2)3/2' 


802 8Di ^ 802 8Di 

dollh ^ 


to get the desired result by ( |3.6[ ), ( |3.7[ ) and ( |3.8[ ). 
'802 


8 t 


80b 

8 t 

802 


+ Ai 


802 

8x 

DiB 


80b 

8x 


802 . 89b 


_^ /rz~^ 

8 D 2 + 8 D 2 ^^ 8 B^ ^ 

802 8 Di , 802 8 Di 

+ 77tV^+Ai 


D 1 D 2 


(1 + 11)12)3/2' 


dDi 8t 


8Di dx 
DiB 




32 


1 /-I I p2 -P 1 -P 2 

1 + 52 ^ ^ (r+|D|2)3/2' 


1 + |D |2 l + lup 

+ ^ ^2 cos ^2 sin 02 ji + 003(02 + 0b) _|_ COS 02 sin 02 P, 

where we used (3.4) for the last two quantities. We now use (3.3) together with the transformation 


( |3.2[ ) to get 

(m _ 80 b \ 

8t 8t ) 
1 


, , (802 80b\ 

^ \ 9a; 9a; / 

: cos^ 02 sin 01 sin 0 B p — 


1 




cos^ 02j 2 — COS 0_B sin 01 sin 02 


+ 


'/l + ^1 

—cos 02 sin02ji + cos(02 + 0 b)- -^ cos 02 sin02p 




l + Di 
iAi 


cos^ 02 sin 9b p — 


1 


i + D'f -“-r- 

D 


l + Di 
cos^ 0232 — COS0B 
D 


— , cos 02 sin 02 — , COS 02p 


+ ^ ^2 cos 02 sin02ji + cos(02 + 0 b) _|_ j 2 COS 02 sin 02P 


pDi 


COS^ 02(cOS02Sin0B + sin02COS0B)- , COS^ 02 

1 + ^? x/ihV^ 

+ cos02 sin02 + cos(02 + 0 b) ^^ 1^2 cos 02 sin02 


l + Dl 

3iDi . 

cos 02 sm 02 — 


1 + i^? 


32 


l + Dl 




COS^ 02 
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<^ 0 ®^^ 2 sin( 6 i 2 + Ob) + cos( 6»2 + 

“r -L/^^ 1 “T -L/j^ 

„• n 32 __„2/I 


1 + i?? 

jiDi a • a 32 

= -^ cos U 2 sin U 2 - , cos t )2 

i+Di 

- y^^^‘2. cos02 (cos 6*2 sin(02 + Ob) - cos(02 + Ob) sin02) 

3l^l „„„ fl a 32 _ ^„2 n pDl _ „ „■ / 


COS 02 sin0B, 


' 1 

ji^i a ■ a 32 2 a 

= rrsfTir^'’’-TTB? 

where we used trigonometric identities several times, and the first equation of (3.51 is obtained. 
dta — {cos(3)da:a = cos 02 ^ ^^^2 si ’2 02 -cos 02 — sin 0 B 


\ ± 

The second equation of (3.5) is obtained by the same 
vector I 2 to (3.1) to get 

/ d02 dD2 dOBdB\ ^ / d02 dD2 BOb dB 
902 DiB 902 . 90b 


- ■ r COS C7‘ 

Vi+m ■ . / 

calculation. Multiply the transformed eigen- 


^9i:>2 dt _, 

902 DiB 902 . 90b f—BB D 1 D 2 

dD 2 + dD 2 ^^ dB^ + (l + |i5|2)3/2/ 


and add the following quantity on both sides: 

Bn. 


902 9L'i 902 dDi _ 902 . 902 


to get 

fd02 , 90b\ , , /902 , 90b\ 

902 DiB 


(JU2 JJilJ 

~~^2 v/TTWVi + B^ 

902 . , , 902 

9^1^! + A 2 9^iP 

DiB 


902 . OOb / 2 D 1 D 2 

(1 + |L»|2)3/2 


p 


+ D,B yrr:^. i ^i^?2 

1 + |iA|2 i + + ^ (i + ii?|2)3/2 

+ --^ cos 02 sin 02 ji - cos (02 - 0 b)^ cos 02 sin 02 p 

1 -\- 33^ 1 3J\ 

1 1 . „ . „ . „ 1 


2 1 2 1 

COS 02 sin 01 sin 0 bP - , cos O 232 + cos 0 b sin 0i sin 62 . 

v/IT^ y/l + I^P 


^0 


l + Dl 


l/l + ^1 V - ■ - 1 

+ . . Lo cos 02 sin 02 ji - cos (02 - 0b) ^ cos 02 sin 02 p 

1 -f- 

: COS 6*2 sin 6*2 — , cos 62 P 


+ 


, ^n 2 ^2 sin ObP - , ^ „ cos^ 02 j 2 + cos Ob , „ ' 

1 + Df \/rD^ \/r+^ 

^ cos 02 sin02ji — cos(02 — 0 b)COS 02 sin 02p 

1 “h 


1 + Dj 


cos^ 02(cos02 sin0B — sin02 cos0b)- , 

1 +^? x/IhL^ 


I 3iDi „ . „ 

+ COS 02 sm 02 

1 + £>1 

jl£>l 


- cos (02 - 0b) cos 02 sin 02 


COS^ 02 

1 


3i^i a • a 32 

TTSf“=''“™''^-yTT^ 


32 2 n 

cos O 2 
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+ 


Dl 


cos 9o sin 9-2 


*^0®^ ^2 sin(6'2 - 9b) - cos(6»2 - 9b) 

1 “T 1 

Jl^l a • a -^2 2 a 

-^ COS 92 Sin 92 - , cos 92 

+ cos 02 (cos 6*2 sin(02 - 9b) - cos(02 - 0 b) sin 02^ 

h _2 n 


jiDi 
1 + Dl 


cos 02 sin 02 — 


./T+m 


cos 02 — 


1 + D( 


J cos 02 


sin( 


Hence, we obtain the second equation of (|3.5|) as follows: 
9t/3 + (cos a)dx(i = cos 02 


jlDl . j2 

- -—77 Sin 62 - , 

1 + Dl a/I + £>1 


0 ■ 0 

cos 02 — Sin 0B 

1 + 


Note that the inhomogeneous term is the same with the first equation of (3.5), and this completes 
the proof. □ 


To summarize, the hyperbolic system (2.6|-(2.7| was rewritten as a quasilinear form (3.1), and 
we transformed it to a diagonal form (3.5) by using the transformation (3.2). 

3.2. Quasilinear hyperbolic systems. Quasilinear hyperbolic systems, especially those which can 
be written as hyperbolic conservation laws, have been extensively studied for last several decades. 
We refer to some good textbooks I1II1I13117] and the references contained therein. In this paper, 
we will use a classical result of Hartman-Wintner ng where local-in-time existence of solutions 
to inhomogeneous quasilinear hyperbolic systems was proved. 

Proposition 3.1. [101 US] t > 0, cc S M, and u = u{t,x) G K". Consider the Cauehy problem 
for the following inhomogeneous quasilinear hyperbolic system with initial data u™: 


(3.9) 

on the following {t^x,u)-region: 

(3.10) 

(3.11) 


Ut = Fux + G, u(0, x) = zi*"(x) 


0 < t < t', |x| < b, 


\ui\ < c, • • ■ 


< c. 


Suppose that the system ( |3.9| ) is strictly hyperbolic on (3.10)-(3.11), where F = F{t,x,u) is a 
matrix function and G{t,x,u) is a continuous vector function on (3.10)-(3.111, and G is of class 
with respect to x, Mi, • • • ,u„. Then there exists a positive t* < t' such that the system (3.9) has a 
unique solution on the region 

(3.12) 0 < t < t* and \x\ -\- Kt < 6, 


where K is chosen such that \\i\ < K, i = 1^ - ■ ■ ,n, on (3.10). 
Proof. We refer to [T6| for the proof. 


□ 


By applying Proposition 3.1 to our case (3.5), we obtain the following corollary. A lower bound 
t* can be determined by following the proof of US); moreover, determining it is much easier than 
that of [ig because the 2x2 matrix in (3.5) is already diagonal. 

Corollary 3.1. Consider the following inhomogeneous quasilinear hyperbolic system (|3.5|).' 


d. 


— cos /3 0 

0 cos a 


cos02(fco sin0B -I- fci sin02 -I- k 2 cos02) 
cos02(fco sin0B -I- fci sin02 -I- k 2 cos02) 


on the following region: 
(3.13) 


|02(t,x)| -b |0B(t,a;)| < 2> 
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where a = 02 — 0b o.nd P = 02 + 0b- Suppose that and /3®” are of class satisfying (3.13) at 
t = 0 and compactly supported in {cc € K : \x\ < P}. Let ki = kpt, x), i = 0,1, 2, he continuous and 
of class with respect to x variable, and assume that there exists a positive nondecreasing function 
k such that 

||fco(0l|oo + ||A:i(t)||oo + ||fc2(t)||oo < Ht). 

Then there exists a unigue solution on a time interval [0,i*) where is defined as follows: 

ft 


t* sup < t : || 02 "lloo + ll^^lloo + 


fc(T) dr < - 


Proof. We note that the system (3.5) is strictly hyperbolic on the region (3.13): 

A 2 — Ai = cos a + cos/3 = 2 cos 02 cosOb > 0. 

The other conditions in Proposition 3.1 hold, hence we have local-in-time existence of solutions 
by applying the proposition. 

In contrast with homogeneous cases, the solution u of (3.9) will grow due to the inhomogeneous 
term as time evolves. For this reason, the domain (3.10) should be restricted as (3.12), i.e., the 
solution u satisfies (3.11) only on the restricted region (3.12), and thus we have to find a time 
interval where the solutions a and /3 satisfy (3.13). We use the integral formula (5.5)-(5.6) in 
Section 5 and estimate them as follows. We first note that 


\02{t,x)\ -f \0B{t,x)\ < - 
By the integral formula (|5.5|), we have 


and \P{t,x)\<^. 


|Q:(t,a:)| < ||a“||oo + / ||fco('r)||oo + ||fci(E)||oo + ||^2(r)||oo 


9 in 11 
S Noo 


-'O 

<ii^^riioo + ii 0 ^ 

Similarly, we have the same estimate for /3. 

m,x)\<\\ 0 r\\oo + \K\\o.+ 


||fco(E)||oo + ||fcl(T)||oo + ||fe(T)||oo dr. 


||*o(e)||oo + ||fcl(E)||oo + ||fe(E)||oo dr. 


Therefore, ( |3.13| ) holds on the time interval [0,3*) by the definition of 3*, and the solutions exist on 
this time interval. □ 


4. Main result 

In the previous section, we derived the system of equations of our interests. The system of 


(2.5)-(2.9) is transformed by (3.2) as follows: 


(4.1) 

(4.2) 

(4.3) 

where a = 02 — ' 

(4.4) 


dtf + vidj + {Di + V 2 B, D 2 - viB) ■ V,/ = 0, 
dta — (cos P)dxa = cos 02(^o sin0B -|- fci sin 02 -I- ^2 cos 02), 
dtP + {cos a)dxP = cos02(fco sin0B -|- ki sin 02 -I- /c 2 cos 02), 
and P = 02 + 0B, and ki, / = 0,1, 2, are given by 


fcn = - 


pDi 

l + DY 


fci = 


jiDi 

l + Dl 


k 2 = - - 


J2 


They are coupled to each other as follows: 

p{t,x) = / f{t,x,v)dv — n{x) and Di{t,x)= / p{t,y)dy, 

J —00 

j{t,x)= / vf{t,x,v)dv, 

JK 2 


(4.5) 
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and the electromagnetic fields D 2 and B are given by 

(4.6) D 2 = + D\ tan 82 and B = tan^s- 

Let /“, 9™, and 0^ be initial data for the unknowns /, 02, and 9b, or equivalently /, a, and /3, 


and we consider the Cauchy problem for (4.1|-(4.6). 


Assumptions on initial data. To study the Cauchy problem for (4.11-(4.6), suitable assumptions 
on initial data should be taken. We will assume that initial data /“, 0™, and 0^ satisfy the following 
conditions A1-3: 

Al. They are compactly supported functions: 

€el{RxR^), 0“eCi(M), 0)j”eCi(K), neCi(M), 

and P denotes the size of their support, i.e., they are supported in 


{(a;,?;) 


|u| < P} or {x : \x\ < P} . 


A2. The fields 0“ and 0^ are bounded such that 


A3. The fields are small in the sense that 


\0 b{^)\ < 2 


for \x\ < P. 


11 ^: 


+ I|6'b Iloo < arctan —. 


Note that there are two kinds of characteristic curves, one from the Vlasov equation (4.1) and the 


other ones from the field equations (4.2)-(4.3). The particles described by the Vlasov equation have 


velocities i), while the Born-Infeld field propagates with speed — cos/3 or cos a, which may be less 
than the speed of light c = 1. Hence, there may be a resonance between a particle trajectory and field 
propagation, and in this case we cannot apply the arguments of m- for the Vlasov-Maxwell case, 
the Vlasov and the Maxwell characteristics are always linearly independent, so x or f derivatives can 
be decomposed into linear combinations of the Vlasov and the Maxwell characteristics (see [131 [E] 
for details). This argument can be applied to our case (4.1 )-(4.3) when we take the third assumption 


above. The assumption A3 implies that the Vlasov and the Born-Infeld characteristics are separated 
at t = 0. Note that 


+ II^B Iloo < arctan ^ 


Qtn 11 
"2 Noo 


+ I|0bIIoo), 


||02”lloo + l|0f?lloo < arctan P < cos 1 

and this implies that under the assumptions Al and A2 we have 

|Di| < cos (max {|q!*"'(x)|, |/3*"(x)|}) = min{cosa*"(x),cos/3™(x)} 

for any x and v satisfying |x| < P and \v\ < P. Remind that 0 “ = 9p — 9“^ and /3“ = 9p+9p^ and 
the cosine function is decreasing on [0,7r/2]. Hence, the particle velocity and the field propagation 
speed are initially separated as follows: 

(4.7) — cos/3“(x) < -Oi < cosQ!“(x). 

We now state the main theorem of the present paper. 


nin 
^2 ; 


Theorem 4.1. Consider the Cauchy problem for the system (4.1 )-(4.6) When initial data f 
and 0g satisfy the assumptions Al-3, then the system (4.1)-(4.6) has a unique solutions f, 02 , 
and 9 b locally in time. 

Proof. The uniqueness is proved by following standard arguments for the Vlasov equations, and we 
only prove the existence part of this main theorem. The following sections are devoted to the proof 
of existence. We refer to Section 5 and 6 for the proof. □ 
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5. A PRIORI ESTIMATES 


In this section, we obtain several a priori estimates for /, a, /3, and their derivatives. The 
arguments basically follow the well-known results [ISIIS]- It turns out that the nonlinearity of the 
field, i.e., the necessity of controlling the field characteristics, complicates the problem and requires 
some restrictions, which will be seen in the proofs of lemmas. 

To obtain a priori estimates for the system (4.1 )-(4.6), we consider a slightly modified version of 
it in this section: for given D* and i?*, we consider the following Vlasov equation: 

(5.1) dtf + vAf + [Dl + V2B*,D* - viB*) • V,/ = 0 


together with (4.2|-(4.6|. We first consider the characteristic equations of the Vlasov and the field 
equations. Suppose that D* and B* fields are on a time interval I. Then, the following system 
of ODEs has a unique solution X{s\t,x,v) and V{s;t,x,v) which are on / x / x K x 


(5.2) 


-X{s) = Vi{s), X{t) = x, 

^V(s) = D*{s,X{s)) + (Ms),-V,{s))B*{s,X{s)), V{t) = v, 


where V(s) = X{s;t,x,v) and V(s) = V{s;t,x,v). The Vlasov equation (5.1) is then solved as 
follows: 

(5.3) f{t,x,v) = r{X{0),V{0)). 

It is well known that the following map is measure preserving: 

{x,v) i-A {X{s;t,x,v),V{s;t,x,v)), 

which implies that ||/(t)||p = ||/™||p for any 1 < p < oo and t G I, and we refer to [H] for details. 

Similarly, if a and /3 are on I, then the following two ODEs have solutions ^(r;t, x) and 
r]{T; t, x) on / X / X M: 

(5.4) :^?(t) = -cos/3(t,^(t)), ^{t) = x, 


(5.5) 


dr 

d 

dr 


? 7 (r) = cos a(r, r]{t) = x. 


Therefore, a and P in (4.2)-(4.3) satisfy the following integral formulae on J x 

pt 


(5.6) 


(5.7) 


a{t,x) = a“(^(0)) 


P{t,x) = nr,{0)) 


/o 


(fco (r, ^(r )) sin 6 »b (r, 5 (r)) -t /ci (r, ^(r)) sin 02 (r, ^ (r)) 
+ (t, ^(r)) cos 02 (t, ^(r))) cos 02 (r, C(t)) dr, 

(^o(r,?7(T))sin0B(T,7?(T)) -t A:i(r, ?7(r)) sin02 (T, 7?(t)) 
-I- fc2(r, 77(r)) cos02(r, 77(r))^ cos 02 (r, 77(r)) dr. 


5.1. A priori estimates for the unknowns and momentum support. We first define the 
momentum support P{t) of / as follows: 

P{t) := sup {|u| : f{s,x,v) 0, 0 < s <t, xGM, v G K^} . 

Note that P{0) = P. For functions D* and B* which are defined on a time interval I, the 
momentum support P{t) is well defined on / because the characteristic equation (5.2) gives 

(5.8) sup \V{t;0,x,v)\ < P + [ ||D*(s)||oo + ||S*(s)||oo ds, tGl. 

Jo 
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The time interval / will be determined later, and we present some useful lemmas. The following two 
lemmas are easily obtained by using (4.4), (4.5), and (5.3). 

Lemma 5.1. Suppose that D* and B* are functions on a time interval I, and consider the 
equations (5.1) and (4.5). Then f, p, j, and Di are well defined and on I and satisfy the 
following estimates on / x M x or / x K.' 


(*) f{t,x,v) <\\P‘ 

{a) \p{t,x)\ < 7r||/“||ooP^(t) + llnhoo- 

(Hi) \j{t,x)\<TT\\f^\\^P^{t). 

(iv) |Di(t,a:)| < ||/“||i + ||n||i. 


Proof. We obtain the first estimate by using (5.3). The second estimate is obtained by the definition 
of p. 


\p{t,x)\< ( f{t,x,v)dv+ \n{x)\ 
JR 2 


<iir 


dv - 


'\v\<P{t) 


<A\r\up\t) + 


The third estimate is similarly proved. For {iv), we use the measure preserving property of the map 
{x,v) i-A (X(0),F(0)) as follows: 


\Di{t,x)\ < 


\p{'t,y)\dy 


< 


[ [ f{t,y,v)dvdy+\\n\\i = \\f^\\i + \\n\\i. 

J —oo Jm? 


□ 


Lemma 5.2. Suppose that D* and B* are functions on a time interval I, and consider the 
equation (4.4). Then ki, i = 0,1, 2, are well defined and on I and satisfy the following estimates 
on I X K; 


(z) \ko{t,x)\ < 7r||/™||ooF2(i) + ||n| 

(a) \ki{t,x)\<Tr\\f'^\\^P'^{t). 

{in) \k 2 {t,x)\ < Tr\\f'^\\^P'^{t). 


Proof. In (4.4), we can see that 

\ko{t,x)\<\p{t,x)\, \ki{t,x)\<\ji{t,x)\, |fc 2 (t,a;)| < |j 2 (t,a;)|. 

Hence, the lemma is proved by Lemma 5.1. 


□ 


In the next lemma, we use the integral formulae (5.6)-(5.7) and Lemma 5.2 to estimate D 2 and 
B fields. Note that a = 02 — Ob and P = 82 + 83 - 

Lemma 5.3. Suppose that D* and B* are functions on a time interval I. Then the system 
(4.2)-(4.3) has unique solutions on I D [0,Ti), which is defined by 


Ti := sup < t : \ \ 0 '. 


9 in 11 
2 \\oo 


+ l |0 


B Moo 


+ / Halloo + SttII/ 


in 11 £>2 

00 


P (r) dr < - 
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and D 2 and B satisfy the following estimates on I D [0,Ti) x K; 

\D2it,x)\ < Vl + (||/“||i + ||n||i)2tan f||6»^"||oo + WSWoo + [ ||n||oo + Stt]|/“||oo-P^(t) dr] , 


|B(i,a:)| < tan |^||6»“||oo + II^B Iloo + J ||n||oo + 37r||/“||oo^’^(r) dr j . 

Proof. The existence interval [0,Ti) is given by Lemma 5.2 followed by Corollary 3.1. Since 62 = 
{a + /3)/2 and 9b = {fi — a)/2, we have the following estimates from (5.6)-(5.7) and Lemma 5.2: 

\e 2 {t,x)\ < l^^"(^(0))l + l/?“(d(0))l ^ ||fco(r)|U + ||fci(r)||oo + ||A:2(r)||oo dr 

< ll^riloo + mWoo + r ||fco(r)||oo + ||fci(r)||oo + ||fc2(r)||oodr 

^0 

< Il^riloo + ||di3"||oo+ f ||n||oo+37r||/“|UP"(r)dr, 

Jo 

and similarly we have for 9b 

\9B{t,x)\ < \\9T\\oo + moo+ f Moo +M\r\\ooPHr) dr. 

do 

Hence, the lemma is proved by ( |4.6| ) as follows: 

|dd2(t,a;)| < + Dl{t,x)tan{\e2{t,x)\) 

< yi + (||/»||i + |H|i)2tan (^||d™|U + ||dij"||oo +£ ||n|U +37r||f"|Up2(r)dr) , 

and 

\B{t,x)\ < tan(|d 2 (t,a:^)|) 

< tan (^lldriloo + II^H Iloo Halloo +3^||/“|UP"(r)dr^ . 

□ 


Remark 5.1. By the same argument with the proof of Lemma 5.3, we can see from (5.6)-(5.71 that 
a, j3, 02 , and 9 b are bounded by the same quantity as follows: 


Qin II ill ntn \ 


\n\\oo + 37r||/“||ooT’^(r) dr. 


max{|Q;(t,a;)|, \f3(t,x)\} < M Iloo -r ||i^b Hoo 
Since a = 62 — 9b and P = 92 -\- 9b, the above inequality gives the following estimate: 

\92{t,x)\ + \9B{t,x)\ < || 6 » 2 ''||oo + II^B Iloo + [ ||n||oo + Stt 1111 ('t) dr. 
In other words, on [0,Ti) we have 


max||6»2(t,a;)|, \ 9 B{t,x)\, |a(t,x)|, |/ 3 (t,x)|| < |d2(t,a;)| + \ 9 B{t,x)\ 


TT 

< 2 - 


Lemma 5.4. Suppose that D* and B* are functions on [0, Ti), and consider the following Vlasov 
equation for V: 

dtp + VMP + (ddi + V2B, D2 - viB) ■ = 0, 
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where D and B are given by the system (5.1) and (4.2)-(4.6), with the same initial data /^(O) = /*' 
Then its momentum support P^(l) satisfies the following integral inequality on [OjTi); 

P\t) < P + 2Vl + (||/“||i + ||n||i)2 

cos-1 (^||0“||^ + ||0“||,^+^'||n|U+3T||r"|UP'(r)dr^ ds, 

where P{t) is the momentum support of f in ( |5.1| ). 

Proof. By Lemma 5.1 and 5.3, we have 
\Di(t,x)\ < ||/“||i + ||n||i, 

\D 2 {t, x)| < \/l + (I I/*" 111 + I |7i||i)^ tan ( I Id™ I loo + 11^^ I loo + [ II' 


+ 3^||/™||ooP"(r)dr , 


\B{t,x)\ < tan |^|| 6 »“ 
Hence, we obtain 

and therefore, 


IW 


B Moo 


+ 37r||/“||ooP"(T)dr . 


|P(i,x)|2< (l + (||/“||i + ||n||if)(l + tan2 

= (l + (ll/“lli + IH|i)^)cos-2(---), 


Ili^Wlloo + ||P(t)||oo < 2^1 + (ll/“lli + IH|i)2 cos-1 
We now use the inequality (15.81) for p and P'^ to get the desired result. 


□ 


5.2. A priori estimates for derivatives of the solutions. In this part, we estimate the deriva¬ 
tives of /, a, and (3. In Lemma 5.3, we introduced a finite interval [0, Ti) on which we could estimate 
the field quantities. To estimate their derivatives, we need an additional interval [ 0 ,r 2 ) which is 
defined as follows: 

T 2 := supjt : lld^lloo + ll^'s Iloo + J ||n||oo+ 37 r||/™||ooP^(s)ds < arctan 



Note that T 2 > 0 because of the assumption on initial data A3. We remind that the assumption 
A3 gives the separation between the Vlasov and the Born-lnfeld characteristics at t = 0. Hence, 
we can see that the interval [ 0 ,T 2 ) is the maximal interval on which the Vlasov and the Born- 
lnfeld characteristics are separated: since a and j3 are bounded by the LHS of the inequality in the 
definition of T 2 (see Remark 5.1), we have on [ 0 ,r 2 ) 


max{|a(t,ai)|, |/?(t,a;)|} < arctan 


1 

wr 


which is equivalent to 


P[t) < cos(max{|a(t,a;)|, |/3(t,x)|}) = minjcosa(t, t), cos/3(t,a:)}. 


Hence, we have 

— cos [3(t, x) <vi < cos a{t, x) 
for any x and v satisfying /(t, x, v) 7 ^ 0 . 

We now fix a closed and bounded interval [0, T] C [0, Ti) n [0, T 2 ), and all the following estimates 
in this subsection will be considered on [0,T]. Note that P{T) is finite, and Ct will denote a 
positive constant depending only on T and P(T) and vary from line to line. We first estimate x and 
V derivatives of /. 
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Lemma 5.5. Suppose that D* and B* areC^ functions on [0,T], and consider the following Vlasov 
equation for p: 

dtp + vAP + {Di + V 2 B, D 2 - nB) ■ = 0, 


where D and B are given hy the system (5.1) and ( |4.2[ )-(4.6 ), with the same initial data p{0) = f™. 
Then we have the following estimate for Vx,vP ■ 

||V,,„/t(t)||oo < ||V,,„/™||oo + Ct [ (l + P,a(s)||oo + ||5,/3(s)||oo)||V,,,/^(s)||oods, 

Jo 

where Ct is a positive constant depending on T, P{T), and initial data. 

Proof. By direct differentiation with respect to x and v, we obtain 
dt{dxP) + vidxidxP) + {D + {v2,-vi)B) ■ \ty{dxp) 

= -{dxD + (v 2 , -vi)dxB) ■ Vyp, 
dt{dy,P) + Vldxidy.p) + {D + {V2, -Vl)B) ■ \t^{dy.p) 

= -idyiVi){dxP) - {dy^V 2 ,-dy.Vl)B ■ Wyp, 


where i = 1, 2, and by taking the characteristic curve (5.2), we have 

rt 


\\^x,vP{t)\\oo < ||V,,,/“||oo+Ct / (l+P,14(s)||oo + PxS(s)||oo)||V,,„/t(s)|Uds, 


/o 


where we used |i?| < Ct by Lemma 5.3. By ( |4.5| ) and Lemma 5.1, we obtain 
(5.9) \dxDi(t,x)\ < \p{t,x)\ < Ct- 


We take x derivative on (4.6) to have 


dxD 2 = 


DidxDi 


tan 62 ~\~ \J 1 ~\~ D‘^ 


dx 02 


pi +Dl '' " COS" U2 

Since \ 92 {t,x)\ < | on [0,Ti), we obtain 

(5.10) \dxD2{t,x)\ < Ct{ 1 + \dxd2{t,x)\) < C't(1 + \\dxa{t)\\^ + ||5a;/3(t)||oo)- 

Similarly, we obtain 

(5.11) \dxB{t,x)\ < CT{\\dxa{t)\\^ + \\dx/ 3 {t)\\^), 

and this completes the proof. □ 


We now estimate the x derivatives of a and j3. In the next lemma, we introduce new quantities 
in order to estimate dxCt and dxd as in [20], but the lemma is proved only when the Vlasov and 
the Born-Infeld characteristics are well separated. We will see in the proof that it is necessary to 
introduce the interval [ 0 ,r 2 ) and restrict all the arguments to it. 

Lemma 5.6. Suppose that D* and B* are functions on [0,T], and consider the following quan¬ 
tities u and w: 

u :={—cos a — cos l3)dxCt and w := {cosacos l3)dx/3. 

Then, they are bounded on [0,T] by a positive constant Ct which depends only on T, P{T), and 
initial data. 


Proof. By direct calculations, we have 

dtu = (sina9ta + suv(idtfi)dxOi + (— cos a — cos fi)dxdta, 
dxU = (sinacla;a + sin f}dxlJ)dxOi + (— cos a — cos/3)i9^a, 
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(5.12) 

dfdxOL + sin/39a;/39a; a — cos/39^a 

= {dxkf) sin^B + fco cosObOxOb + dxki sin02 + cos025a;^2 + dxk 2 cos02 — ^2 sin029^02) cos 02 

— (fco sin0B + fci sin02 + /c2 cos02) sin029^02- 

We now derive an equation for u as follows: 


dtu — cos jSdxU 

= sinadtadxOi + sin/39(/39a;a — cos adxdta — cospdxdta 

— cos /3 sin a{dxa)‘^ — cos /3 sin (3dxPdxCt + cos a cos /39^a + cos^ /39|a 
= suYadxa{dta — cos jSdxOi) + s\nj3dxa{dtj3 + cosa9a/3) 

— (cosa + cos fi){dtdxOi + sinjidxfidxa — cos/39^a). 


Since the inhomogeneous terms of (|4.2[) and (4.3) are same, we have 


dtU — cos PdxU 

= (sin a + sin l3)dxCt{ko sin0 b + ki sin 02 + fc 2 cos 02 ) cos 02 

— (cosa + cos/3)(9a;fco sin0B + dxki sin02 + dxk2 cos02) cos02 

— (cosa + COS/3)(fco COsObOxOb + ki COS029a;02 — k 2 sin029a;02) cos 02 
+ (cosa + cos/3)(A:o sin0B + fci sin02 + fc 2 cos02) sin029a;02 

=: + ^2 + /a + I 4 , 


where we used (5.12). For /i, we have 


^ sin a + sin /3 ,, . ^ ^ ^ , 

ii =- 7 ;u{ko sm 0 B + fci sm 02 + fe 2 cos 02 ) cos 02 

cos a + cos p 

= —u{kosm 02 sin 0 B + ki sin^ 02 + ^2 cos 02 sin 02 ) 

= — ^^fco(cosa — cos/3) + fci(l — cos( 202 )) + A :2 sin( 202 )^, 


where we used trigonometric identities. For / 2 , we have 

I 2 = —2 cos 02 COS0B(9a;fco sin0B + 9a;fci sin 02 + dxk2 cos 02 ) cos 02 

= — 2 (C 0 S^ 02 COS 0 B sin 0 B 9 a;fco + COS^ 02 sin 02 COS0B9a;fci + COS^ 02 COS0B9a;fc2)- 


(5.13) 


For J 3 and I 4 , we note that 

M + W 


- li 

2 = (cos a + cos j3)dx9B and —^— = (cos a + cos /3)dxd2- 


Hence we have 


r f, „ u + w , ^ w — u , . ^ w — u 

h = - [ko cos Ob — h ki cos 02 — -fc 2 sm 02 —- — 

XV — Xi 

I4 = —-—(fcg sin 0 B + ki sin 02 + k2 cos 02 ) sin 02 , 


cos 02 


and then J 3 + I 4 is written as follows: 


h + h 

'll / \ 'uj / 

= - - fco cos a + fci cos(202) - ^2 sin( 202 ) j + 


fcg cos /3 


fci cos( 202 ) + /c 2 sin( 202 )^ . 
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We write Ii+ I ^-\-14 as follow for simplicity: 

+ ^3 + h 

1 


= u fcp — cos a + - cosP I + ^1 + cos(202) I — ^2 sin( 202 ) 


+ w ^-i/coCOS^- cos(26i2) + ^fc2Sill(202) 

=: Kiu + K2W, 

where Ki and K 2 are elementary functions of fco, fci, ^ 2 , a, and f3. We now control the derivative 
terms dxki, i = 0,1,2, in 12 - By direct calculations, we have 


(5.14) 

and similarly, 

(5.15) 

(5.16) 


Oxkn = - w 

^ ° 1 + Df 


Di 

l + Zl? 


dxP- P 


^xP P 


1 

(1+^ 

1 - Dl 

(T+W’ 


{dxD4{l + Dl)-2DldxD4) 


o , _ o • , ■ 1 - 


dxko — 


y/IT^ 


dxj 2 + j 2 P 


D, 




If we plug the above results into (5.13), then we can see that I 2 is written as follows: 

I2 = Vt^dxP + ^idxji + ^2dxj2 + Lq Li + L2, 

where fli, i = 0,1, 2, are elementary functions of a, /?, and Di, while Li, i = 0,1, 2, are elementary 
functions of a, /3, p, ji, j 2 , and Di. Consequently, we obtain the following equation for u: 


(5.17) 


dtu — cos pdxU = Kiu + K 2 W + LlodxP + i^idxji + Ll 2 dxj 2 + Lq Li ^ 2 - 


By Lemma 5.3, the system (4.2 )-( |4)3| ) has unique solutions on [0,T], hence the characteristic 
curve (5.4)-(5.5) are well defined on [0,T]. Along the curve (5.4), we can rewrite (5.14) as follows: 

ct 


(5.18) 


w(t,a;) = u“(^(0)) + f (Kiu + K 2 W + Lq + Li + L 2 )iT,^{T)) dr 
Jo 

+ / {LlodxP + Llidxji + LL 2 dxj 2 ){T,i{T)) dr. 

Jo 


For simplicity, the last three term of (5.15) will be denoted by 


J, = / n,{T,^{T)){dxji)iT,^{T))dT 

Jo 


for i = 0,1, 2 and jo = p. 


To estimate each Ji, we use Glassey-Strauss’ argument in m- We first introduce the differential 
operators T and S: 

T := dt — cosjJdx and S ■.= dt + vidx, 
by which x derivatives can be written as 

S-T 
vi + cos /3 ’ 


a 



















THE VLASOV-BORN-INFELD SYSTEM 


21 


Note that the denominator above does not vanish on [0,T] because we have — cos/3(t,a:) < i)i < 
cos a(t, x) on [0, T 2 ) for x and v satisfying f{t, x, v) 7 ^ 0. In other words, the assumption A3 and the 
restriction to [ 0 , 72 ) enable us to apply the argument of [TH]- Jo term is estimated as follows: 

Jo= [ f2o('r,^(T))(a^p)(r,^(T))dr 
Jo 

= [ [ ^o{T,^{T)){dxf){T,^{T),v)dvdT - f no{T,^{T)){da;n){^{T))dT 
Jo 7 r 2 Jo 

= '■ Joi + >^02- 

We rewrite Jqi in terms of T and S as follows: 


>^01 — 


S-T 
Di + cos (3 
1 


[ [ f^o('r,C(T)) 

Jo^ 7k2 

Jo +cos/3(T,^(r)) 

- f f ^o{t,^{t))- 

Jo 7r2 1 


0 jR2 wi+ cos/3{t,^(t)) 

^^o(T,^(r)) ^ 


fiTAiT),v) dv dr 

{Sf){T,^{T),v) dv dr 
(Tf)iT,^{T),v)dvdT 


0 7r 2 hi +COS/3 (T,^(r)) 

^o('r,C(T)) 


(5.19) 

/ / V, 

(5.20) 

Jo 7 r 2 

(5.21) 

-/ 


7r2 \v 

We remind that 



0 7 r 2 hi +cos^(t,^(t)) 
^^o(r,C(r)) 


dr 


(-77 - (h 2 , -vi)B)f{T, ^(t), v) 


f{T,^{T),v) 


dv dr 


a 


hi +cos/3(t,^(t)) 


hi +cos/3(T,^(r)) 


dv dr 

(77 + (h 2 , -hi)B)/(r, ^(r), v) dv dr 
fi'r,^iT),v) dv dr 


;fit,x,v)-^, ) dv. 


hi + cos/3(t, a;) ’ ’ hi + cos/3“(^(0))' 

2 Di 

f 2 o = 2 cos 02 cos 0 _B sin^Bi- 

1 + 77f 

and we need the followings, which can be verified by direct calculations: 

a, 


1 

1 1 

(5 ^ 

hi + cos/3 

+ [up (hi +cos/3)2 

1 1 + klV 


dr 


6'2('r,C(T))J = --w(t,^(t)). 


Sb{t,C{t)) =+inhom, 


dr 


DiiT,C{T)) = -ji(r,C('r)) -cos/3(t,C(t))p(t,C(t)), 


where ‘inhom’ is the inhomogeneous term of (4.2) evaluated at (t,^(t)). Note that 77, B, p, j, ko, 
ki, and k 2 are bounded by a constant Ct on [0, T] by Lemma 5.1-5.3. By applying the above results 
to (5.16)-(5.18), we obtain the estimate for Jqi. 

iJoil < Ct [ [ -rr—- 1, ,,,, dvdr 

Jo J\v\<Pir) ki+cos/3(r,C(T))|2 

B r b/iio" 

dv dr 


Ct 
C , 


/o 7kl<P(r) 1^1 +cos/3(t,^(t))|2 
^ —ATT—iT dp C 

\v\<p{t} 1^1 +C0S/7(t,a:)| 
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^ Ct H“ Ct 


[ IkDIIc 

Jo 


, dr. 


In the last inequality, we used the fact that — cos/?(s,x) < Di < cosa(s,a;) on [OjTa) for any x and 
V satisfying f{s,x,v) ^ 0. J 02 is clearly bounded by Ct'- \Jo2 \ < Ct, and Ji and J 2 are estimated 
by the same way: 

rt 


\Ji\ < Ct + Ct 


[ lk(7 
Jo 


, dr, 


z = 1,2. 


The other terms, Ki, K 2 , To, Ti, and L 2 in ( 5.14| ) are easily bounded by Ct on [0,T]. Therefore, 
we obtain the following estimate for u from (5.151: 

\\u{t)\\ao < Ct + Ct [ ||u(r)||oo + ||w(r)||oo 

Jo 

In the same way, we have the following estimate for w after long calculations: 

\\w(t)\\^ < Ct + Ct [ ||m(t)||oo + ||w(t)||oo dr. 

Jo 

We apply the Gronwall inequality to obtain the desired result, and this completes the proof. □ 


Lemma 5.7. Suppose that D* and B* are functions on [0,T], and consider the following quan¬ 
tities: 

||9xa(t)||oo, 11 I loo, \\dxD{t)\\^, \\dxB{t)\\^, ||Va;,„/'l'(t)||oo, 

where f^ is given by Lemma 5.5. Then, they are bounded on [0,T] by a positive constant Ct which 
depends only on T, P(T), and initial data. 

Proof. By Lemma 5.6, we have the following estimate on [0,T]: 

Ct > {cosa{t,x) + cos I3{t, x))\dxC({t, x)\ = 2cos92{t,x) cos0B{t,x)\dxa{t,x)\. 


On the other hand, we have on [0,ri) (see Remark 5.1) 


\& 2 it,x)\<^, \9B{t,x)\<^. 


Therefore, we obtain on [0, T] 


cos92{t,x) > cos9Bit,x) > 

Oj’ Oj’ 

and this gives ||92,Q!(t)||oo < Ct- Similarly, we obtain the boundedness of ||cla;,d(i)||oo- The bound¬ 
edness of dxD and dxB is obtained by (5.9)-(5.11) as follows: 

||5,i?(i)||oo + ||9.S(t)||oo < Ct( 1 + ||9.a(t)||oo + P,/3(t)||oo) < Ct- 


By using Grdnwall’s inequality on the result of Lemma 5.5, we can see that x,vfH'^)\\oo is bounded 
by Ct- This completes the proof. □ 


6 . Iteration scheme and proof of convergence 


In this section, we prove the main theorem. We first introduce the iteration scheme for (4.1 )-(4.6) 
and show that the sequence converges to a function by using the a priori estimates obtained in 
the previous section. 
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6.1. Iteration. We use a standard iteration. We take sequences /”, 6*27 0^, n = 1,2, • • •, as 

follows. Define f := /“, 0^ := 6»^”, and 0% := 0™, or equivalently := - 0“ and 

/3° := /3*” = 0“ + given (n — l)-th step, we define n-th step as follows: /" is taken as the 

solution of the following equation: 


( 6 . 1 ) 


dtr+vAr + iD^-^+V2B 


,D 


n —1 


= 0, 


rio,x,v) = r{x,v), 

from which we take p", j”, and D" as follows: 


( 6 . 2 ) 


p"(t,a;)= f f'^(t,x,v) dv — n{x) and Di{t,x) = f p^{t,y)dy, 

J—00 

j^{t,x)= [ vP{t,x,v)dv, 

JR2 


and then we get /c”, i = 0,1,2, from p”, j”, and D” by using (4.4). 02 and 0g are taken as the 
solution of the following inhomogeneous quasilinear hyperbolic system: 

i9ta" — (cos/?")9a;a" = cos 02 (feg sin0^ + fc" sin 02 + fcj cos 02 ), 

(6.3) i9t/3” + (cosa”)9a;/3" = cos02 (feg sin0)^ + fc" sin02 + fcj cosd^), 

02 (0, x) = 02 "(a;) and 0^(0, a;) = 0^(a;), 

where a" = 02 — 0^ and /3” = 02+0)^. Finally, and B" are obtained: 

(6.4) D” = ^1 + (D”)2tan0J and B”=tan0g, 
and this completes the n-th step. 

6.2. Some remarks. We remark that the a priori estimates obtained in the previous section can 
be applied to each n-th step uniformly on n. In the previous section, we Hrst defined the momentum 
support and then obtained an integral inequality for the momentum support in Lemma 5.4. We now 
redefine B(t) as the solution of the following integral equation: 

F(t) = F + 2v/l + (||/“||i + ||n||i)2 

(6.5) F ( /•» 

ll^riloo + ||0)3"||oo+ / ||n||oo + 3^||/“||ooP'(r)dr ) ds, 


X / cos 

Jo 


-1 


where P, /*", n, 0“, and 0^ are given initial conditions satisfying the assumptions A1-A3. We 
dehne a time interval [0, To) as the maximal interval on which the solution of (6.5) exists. Then, Ti 
and T 2 are also redehned as follows: 

pt 


Ti := sup t : || 0 riloo + \\ 0 ' 

+ \\o: 


tin 11 

B Noo 


.in 11 
B Noo 


+ 3^||/“|UP"(s)ds< - 


+ 37 r||/™||ooP^(s) ds < arctan 


T 2 := supU : ||02"|| 

I ./o 

We now take 

/:= [0,To)n[0,Ti)n[0,T2). 

On the other hand, the momentum support of is defined as follows: 

P„(t) := sup {|u| : /"■'■^(s,x,u) 7 ^ 0 , 0<s<t, ccGK, v £ . 


F{t) 


If we compare the result of Lemma 5.4 with (6.5) and use the mathematical induction on n, then 


we can see that each P„(t) is bounded by F(t), i.e., 

( 6 . 6 ) P„(t) < F{t) 
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for any n. This implies that the existence interval of P„, say [0,To „), contains that of P, i.e., 

(6.7) [0,To) C [0,To,„). 

Consider now the iteration functions /", o", and /?". The 0-th iteration functions are clearly 
on 1. Assume that n-th iteration functions are of class on I. Then, _D" and i?" are C^, and the 
solutions X” and of the following characteristic system exist on I, 
d 

( 6 . 8 ) 

ds 

and therefore /"‘''^(t, x, v) = /“(X"(0), C"'(0)) is a solution on I. We now obtain (n l)-th step 
of the quasilinear hyperbolic system ( |6.3[ ) with i = 0,1, 2, which are on I and satisfy 

\\K^\moo + ||fcr'(i)||oo + ||fc 2 +'(t)||oo < \\n\\^+3n\\r\\^P^{t). 

Hence, by Corollary 3.1, the (n -|- l)-th step of (6.3) has solutions on a time interval [0,ri^„) 
which is defined by 


W"(s) = Ci"(s), X^(t)=x, 

-l/”(s) = P"(s,X"(s)) + ( 14 "(s),-Ci"(s))i?"(s,X"(s)), C"(t) = n, 


Ti,„:=sup t:||0riloo + ||0j?||oo+ / IU + 37r||/“|UP„"(s) ds < w 


By the definition of Ti and (6.6), we can see that 
(6.9) 


[0,Ti) c [0,Ti, 


and therefore and are solutions on I. P"+^ and are clearly on / due to (6.2) 
and (6.4), and this completes the (n + l)-th step. To summarize, the iteration functions are well 
defined as functions and exist on I uniformly on n. 

We finally consider the separation between the Vlasov and the Born-Infeld characteristics for each 
n-th step. A time interval [0,T2_„) is defined as we did in Section 5.2. 

rt 


T2,n :=supa: ||dri|oo + ||d: 


B Moo 


|n||oo -l-37r||/“||ooP^(s)ds < arctan 


1 


Pn{t) 


By (6.6), we can see that 

( 6 . 10 ) 

and by the same argument we have 

— cos 


[0,T2)C [0,T2,„), 

?"^^(t,a;) < Di < cosa"^^(t,a;) 


for any x and v satisfying /”+^(t,a;,n) yf 0 on [0,T2) for any n. Hence, on [0,T2) the Vlasov and 
the Born-lnfeld characteristics of each n-th step are well separated uniformly on n. 

Consequently, we have (|6.7|) 


_ and ( |6.10[ ), i.e., I C [0,Ti^„) n [0,r2,„) C [0,r3_„) for any n, 

and all the a priori estimates obtained in Section 5 are applied to iteration functions on I. We now 
fix a closed and bounded interval [0,T] C I, and on this time interval the following quantities are 
bounded by a positive constant Ct which does not depend on n: 


( 6 . 11 ) 


V I nn 02 n 02 -r 

X,vf 5 ^xP 5 ^xJ 


dl 


dir, dlD-, dlB^ 


where f = 0,1 and n > 0 are integers. 


6.3. Proof of the convergence. We now show that the iteration functions {/"}, {«"}, and {/!”} 
are Cauchy sequences in C^. The proof is straightforward and will be given by a sequence of 
lemmas. We divide it into two parts: the estimates for the iteration functions themselves and for 
their derivatives. Remind that we fixed a compact interval [0,r] C I, and all the following lemmas 
and arguments will be considered only on [0,T]. The quantities in (6.11) will be roughly estimated 
by Ct- 
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6.3.1. Estimates for the iteration functions. We first show that the sequences {/"}, {«"}, and {/?”} 
converge to continuous functions. 

Lemma 6.1. Consider the iteration functions (6.1)-(6.4|. Then we have 

wnt) - rmoo < Ct f\\D^-\s) - D--\s)U + \\B^-\s) - B^-\s)\U ds. 

Jo 

Proof. For any n and m, we have 

dtr+vid.r + + (^2, • v,r =«, 

dtr + vAr + + {V 2 , = 0 , 

and by direct subtraction and using ( |6.8[ ) we have 
\p{t,x,v) - f"^{t,x,v)\ 

< {\D^-\s,X^-\s)) - D^-\s,X^-\s))\ 


+ |B™-i(s,X"-i(s))-B”-i(s,X”-i(s))|)||V„/'"(s)||oods. 
By taking supremum, we get the desired result. 

Lemma 6.2. Consider the iteration functions (6.1)-(6.4|. Then we have 
\\D^{t) - D^{t)\\^ + ||B™(t) - B^{t)\\^ 

< CT{\\rit) - rmoo + \m{t) - a^moo + wmt) - 

Proof. By (6.2), we have 

DT{t,x) - DT{t,x) = [ [ f"^{t,y,v) - r{t,y,v)dvdy. 

J — oo 


We use (6.6) to have 

( 6 . 12 ) 


mrx) - D^{tm\ < crwrr - rmoo. 


and 


□ 


By (6.4), we have 

\D^rx) - Dtrrx)\ < crwrit) - riru + iKit) - erru, 

\B^it,x) - B^{t,x)\ < CrWO^it) - 0SWIIoo. 

Note that 

(6.13) \\9T{t) - enrioo <\m{t) - armoo + wmt) - rmoo, * e {2,b}, 

and therefore we get the desired result. □ 

Lemma 6.3. Consider the iteration functions (6.1)-(6.4|. Then we have 

2 

j2\\kT{t) - rmioo < Crwrit) - rrm 

Proof. Note that 

\jr{t,x)-r{t,x)\< [ \r{t,y,v)-r{t,y,v)\dv<CT\\r{t)-rmoo, ze{o,i,2}, 

where we used ( |6.6[ ), and p" was denoted by jo for simplicity. Together with (6.12), we obtain the 
desired result. □ 
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Lemma 6 . 4 . Consider the iteration functions (6.1)-(6.4|. Then we have 

\\a^{t)-a^{t)\U + \\r{t)-P'^moo 


< 


Ct f Wa^ir) - a"(r)|U + \\r{r) - /3"(t)|U + ||r(r) - /"^lU dr. 
Jo 


Proof. We use (6.31 for m and n. 

dta'^ - [cos l3'^)d^or = cos 9^ [kff sin Off + /fc™sin6»™ + cos6»^), 

dta'^ - [cos I3^)d^a^ = cos6»”(/c” sin6»g + sinOtf + fcj cosOtf), 
After direct subtraction and applying the following characteristic curve: 


(6.14) 


_^-(r) = -cos/3™(T,r(r)), r[t)=x, 


we make the same argument with Lemma 6.1. We use (6.13) and Lemma 6.3 to have 
\a^[t,x) - a^[t,x)\ 

< Ct f \\or[T) - a"(r)|U + ||/3™(t) - /3"(t)|U + ||r(r) - r(T)|U dr. 

Jo 

By the same calculation, we obtain the estimate for /3™ — /3", and this completes the proof. 


□ 


Convergence to continnons fnnctions. By Lemma 6.1, 6.2, and 6.4, it is proved that the 
sequences {/"}, {o"}, and {/3"} converge to continuous functions. For simplicity, we set 

r^[t) :=\\r[t)-r[t)\u 

9^-[t) := ||a™(t) - a"(t)|U + \\r[t) - /3”(t)||oo. 

Then, the above results are written as follows: 

(6.15) /"*”(<) <Ct [ 

Jo 


(6.16) 


e"^^[t) <Ct f /’""(s) + 9^^[s) ds. 
Jo 


We apply Gronwall’s inequality to (6.16) to have 
(6.17) 


9"^^[t)<CT [ r"[s)ds, 
Jo 


and then apply it to (6.15) to have 



/ 

mn(^) ^ 

Ct 


+ 





lo 

pt 

Jo 



= 

Ct 


ds + Ct 

(6.18) 


< 

Ct 

fo 

ds. 

By iterating it. 

we obtain 






jmn 

^(i)<^(CTT)'=||/“|U : 


for m,n > k, 


where the constant Ct is the same one in the last inequality of (6.18). This implies that {/"} is 
a Cauchy sequence in C° norm, so are {a"} and {/?"} due to (6.17). Consequently, the iteration 
functions converge to continuous functions /, a, and /3. 
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6.3.2. Estimates for their derivatives. In this part, we estimate the derivatives of /", a", and /3" 
to confirm that the solution is C^. We introduce a notation £„„ for a small positive quantity which 
depends on n and m such that £nm tends to zero as n and m go to infinity. This quantity may 
depend on T but not on the other variables, and its value will vary from line to line. We also use a 
handy notation for {x,v) € x as z = {x,v) € Note that 

kl = |(a;,n)| < |a:| + |z;| < Vii\z\, 

where | • | is the usual Euclidean norm on For simplicity, Z" will denote (X”, E") in some places, 
i.e., 

Z'^{s;t,z) = (X"(s; t, T, u), F"(s; t, x, t>)) S K^. 

In this part, we follow the arguments and notations of |18j . 


Lemma 6.5. Consider the characteristic system (6.8). 


-x"(s) = x-(t) = x, 

as 

Ae-(s) = F^is, X"(s), E"(s)), V^it) = V, 

as 

where F^(t, z) = D'^{t, x) + {v 2 , —vi)B^{t, x). For any n ^ m, we have the following estimates: 


(*) \X^{s]t,z) - X"'{s-,t,z)\ + z) - V"^{s;t,z)\ < £nm- 

{ii) |V,X"(s;t,z)| + |V,y’^(s;t,z)| < Ct- 

{ill) I V,X"(s; t, z) - t, z)| + | t, z) - t, z)| 

< Snm + CT [ |(a.E”)(r,X”(T)) - (9,E”)(r,X"^(r))| + - a,F™(r)|U dr, 

J S 

where £mn is a small positive quantity which tends to zero as n, m —> oo and depends on T but not 
on X and v. 


Proof, {i) We consider the characteristic systems for n and m together with the same initial data, 
Z^{t) = Z™‘{t) = z. By integrating the systems from s to t, we have 

rt 


and 


|X"(s;t,z)-X-(s;t,z)|< [ f, z) - E-(r; t, z)| dr, 

J S 

|l/"(s;t,z)-y™(s;t,z)| 

< f |F”(T,X”(T))-E”(T,X’"(r))| + |E"(r,X™(T))-F™(T,X'"(r))|dr 

J S 

< Ct / |X”(r;t,z)-X”"(T;t,z)|dr + enm, 

J S 

where we used the fact that is bounded and {F”} is a Cauchy sequence. We combine the 

above two inequality to have the following estimate: 

|X"(s;<,z)-X’"(s;t,z)| + f |X”{r; t, z) - X™(t; t, z)| dr, 

J S 

where the constant Snm is another small constant depending on T, but it still satisfies the property 
in the statement of the lemma. By applying the Gronwall inequality again, we obtain the desired 
result for X”, 

|X”(s;t,z) - X”"(s;t,z)| < < £„„, 

and for E", 

|E”(s;t, z) - V’^{s;t,z)\ < Enm, 
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and this proves the first result. 


(ii) By direct differentiation with respect to x, we obtain 

(6.19) =Vv^Vr ■ d,V^is;t,z), 

(6.20) a.t>”(s; t, z) = (5,L»”)(s, X"(s; t, z))5,X"(s; t, z) 

+ -Vv. t/i")5,F"(s; t, z)B^{s, X”(s; t, z)) 

+ i%^,-KmB-)is,X-is;t,z))d,X-is;t,z), 

where the dots denote s derivatives. Since dxX'^{t) = 1 and dxV^{t) = 0, z = 1,2, by integrating 
from s to t we obtain 


|9x^”(s)| < l + C* / |9,y"(a)|da, 

J S 

\d.V^{s)\ < Ct f |5,X"(a)| + \d,V^ia)\da, 

J S 


and then Gronwall’s inequality gives the desired result after we estimate v derivative quantities by 
using the same arguments. 


(in) To prove the third estimate, we use (6.19)-(6.20) for n and m. Note that the map u i—> f) is a 
smooth function with a bounded norm. 

\d,x^{t) - d,x^{t)\ 

< f |Vy.l>i"-Vv™l/ri|5.F"(r)| + |Vv^l>ri|5.F”(r)-5,y“(r)|dr 

J S 

<enm + C [ |9,F”(t) - 5,F™(t)| dr, 

J S 

where we used (i) and (ii). By the same way, we have 

<enm+CT f \(d,D^){T,X^(T))-(d^D”^)iT,X"^iT))\ + \d,X^(T)-d,X^(T)\dT 

J s 

+ Ct r |5,F”(r) - a,y™(r)| + |i3"(r, X"(r)) - B^{t,X^(t))\ dr 
J S 

+ Ct r|(9.il")(r,X"(r)) - (a,il™)(r,X-(r))| dr, 

J S 

where we used (i) and (ii) together with the fact that the quantities in (6.111 are bounded. In the 
above inequality, we can see that 

|(a,i^")(r,X"(r)) - (d^D^)(T,X^(T))\ 

< \(d,D-)(T,X-(T)) - (a,i^")(r,X“(r))| + \(d,D-)(T,X^(T)) - (d,D^)(T,X”^(T))\ 

< |(a,l?”)(r,X”(r)) - (d,D^)(T,X^(T))\ + P,D"(t) - d,D^(T)\\^. 

By the same arguments, we have the similar estimate for dxB'^ — dxB^. 

\(dxB-)(T,X^(T)) - (a,i?-)(r,X™(T))| 

< |(5,B")(T,X”(r)) - (dxB^)(T,X^(T))\ + ||9.il"(r) - dxB'^(T)\\^. 





THE VLASOV-BORN-INFELD SYSTEM 


29 


For i?" — term, we have 

|B"(T,X"(r))-B™(r,X-(T))| 

< |B"(t,X”(t)) - B"(r,X™(r))| + |B”(t, X'"(r)) - X™(r))| 

^ ^nm ; 

where we used (i), boundedness of and the Cauchy property of {B"}. We combine the above 

results to obtain 

\da,V^{t) - d,v^{t)\ 

< enm + Ct f |a.W"(T) - dr + |a,F"(r) - d,V^{T)\dr 

J s 

+ Ct f |(a,F")(r,X-(r)) - (5,,B")(r,X™( t))| + ||a,B"(r) - a,B-(r)|Udr. 

J s 

By applying Gronwall’s inequality, we obtain the estimate for x derivatives. The estimate for v 
derivatives is verified by the similar calculations, and we obtain the desired result. □ 


Lemma 6.6. Consider the following two characteristic equations: 

(6.21) |_^-(^) = _cos/3"(r,r(T)), Cit) = x, 

(6.22) ^ cosa”(r,r 7 "(T)), ? 7 ”(t) = x. 

For any n m, we have 

(i) \C{T-,t,x) - £,'^{T;t,x)\ + \r]'^{T;t,x) - 77™(t; t,a;)| < £„„. 

(a) \d^C{x;t,x) \ + \d^r]"-{T-t,x)\ < Ct- 

{in) \d^C{x;t,x) - d^f'^{T;t,x)\ 

< Snm + CT J |(a,/3")(s,r(s)) ” (9.^")(s, C(s))| + \\d,l 3 ^{s) - {s)\\oo ds. 

(iv) \dxiq"'{T;t,x) - d^r]""{T;t,x)\ 

< £nm + Ct f |(9a:a”)(s, ? 7 ”(s)) - (Sa;^”)(s, T?™(s))| + \ \dj;a^{s) - 9a;a™(s)||oo ds. 


where £mn is a small positive quantity which tends to zero as n,m ^ oo and depends on T but not 
on X and v. 


Proof, (i) The proof is almost the same with that of the previous lemma. We first consider the 
characteristic equation (6.21) for n and m together with the same initial data C"(0 = C™(0 = 

By integrating the equation from r to t, we have 


\C{r-,t,x)-r{r;t,x)\ 

< ^ |cos/3"(s,r(s)) - cos/3’-(s,r(s))l + | cos/?™(s,r(s)) - cos/3'”(s,r(s))| ds 
<£nTn+CT J \C{s;t,x) - f'^{s;t,x)\ds, 

where we used the fact that {/?"} is a Cauchy sequence and dx/3‘^ is bounded. By applying the 
Gronwall inequality, we obtain the desired result for 


\C{x-,t,x) - ^”"(r;t,x)| < £nm 
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We apply the same argument to (6.22) to obtain the second estimate for ry”. 


{ii) By direct differentiation with respect to x, we obtain 

(6.23) = sin/3"(T,r(r))(9./3”)(r,r(r))5,r(T), 

which gives 

\d,C{r)\<l + CT j\d,C{s)\ds, 

and we obtain the desired result by Gronwall’s inequality. The estimate for is given by the same 
way, and we skip it. 

(in) The proof of the third estimate is almost the same with that of Lemma 6.5 (in): we use 


(6.23) instead of (6.19)“(6.20), consider the equation (6.23) for n and m, and use the result (i), 


boundedness of dxP'^, and the Cauchy property of {/^”}, and then we obtain the desired result. 

(iv) The last estimate is similarly verified as in (Hi), and we skip the proof. □ 


Lemma 6 . 7 . Consider the iteration functions (6.1)-(6.4). Then we have 


\\dxF^(t) - < enm + - a,/ 3 "(t)||oo), 

where F"^(t,x,v) = D'^(t,x) + (v 2 ,—Vi)B^(t,x) and Smn H a small positive quantity which tends to 
zero as n,m ^ oo and depends on T but not on x and v. 

Proof. We note that the first component of is written as follows: 

dxD-^ - dxD^ = p™ - p" = j r-r dv. 

Since {/"} is Cauchy and the momentum supports of /" are bounded by P(T) uniformly on n, we 
have 

\dxDT -dxD^\<er,m- 


We use (6.4) to estimate the second component of — d^D'^, i.e., 

'a”+/3” 


tan 


+ ^l + (D-)2- 


1 


5a; a” + 5a; 5" 


COS' 




yiTTBrp 

By direct calculations, we have 

\dxD^-dxD^\ 

1 

<CtJ 2 {\91DT - dlD^\ + \d>^ - + \dir - dm) 

mo 

< £nm + C'r(|5a;a'" - 5a,a"| + |5a;/3"‘ - 5a;5"|), 
and the same argument gives the estimate for dxB'^ — dxB'^. 

\dxB^ - dxB^\ < + Cr(|5,a™ - 5,a"| + |5,/3”^ - 5,/3"|). 

This completes the proof of the lemma. 


□ 
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Lemma 6 . 8 . Consider the characteristic system (6.8). 


x'^{t) = x, 

Al-(s) = F^is, X-is),V-is)), V^it) = V, 

where F"‘(t, z) = D'^{t, x) + (02, —vi)B^{t, x). For any z ^ z', we have the following estimates: 

{i) |X"(s; t, z) - X"(s; t, z')\ + \V^{s; t, z) - V^{s; t, z')\ <Ct\z-z'\. 

(ti) \X,X^is-,t,z)-X,X^is;t,z')\ + \X,V^{s;t,z) - V,F"(s;i,z')| 

<Ct\z-z'\+Ct f \id,D^){a,X^ia;t,z)) - id,D^){a,X^ia-t,z'))\da 

J S 

+ Ct f \{d.B^){a, X”(a; t, z)) - {d,B^){<j, X"(a; t, z'))\da, 


where the constant Ct does not depend on n. 

Proof, (i) We apply the mean value theorem to Z^{s; t, z) — Z^{s; t, z') and use Lemma 6.5 (zi), and 
then the first estimate is obtained. 


{ii) Let 2 ^ z', and consider two characteristic curves for {a;t,z) and {a;t,z'). By using (6.19)- 
(6.20), we have 

\dj;X'^{s;t,z) - dxX'^{s;t,z')\ < Ct\z - z!\ + Ct f Id^V^-ia-t, z) - dxW^{a;t, z')\ da, 

J S 

where we used the result (z) and Lemma 6.5 (zz). In the same way, we obtain 
\d,V^{s;t,z)-d,V^{s-,t,z')\ 

<Ct f \{d,D^){a,X^ia;t,z)) - id,D'^){a,X'^ia-,t,z'))\da 

J S 

+ Ct [ \d,X^{s-,t,z)-d,X^{s;t,z')\da 
J S 

+ Ct f \z - z'\ + \dxV'^{s;t,z) - dxV'^{s;t,z)\da 

J S 

+ Ct f mB^){a, X'^l^a-1, z)) - {d,B^){a, t, z'))\da, 

J S 

where we used (z). Lemma 6.5 (zz), and the fact that d^D^, B^, and d^B^^ are uniformly bounded. 
Gronwall’s inequality again gives the desired result after we estimate v derivative quantities by the 
same arguments. □ 

Remark 6.1. We can see that dj;D^{t, x) — dxD”‘{t, x') and dxB'^{t, x) — dxB'^{t, x') are estimated 
by the transformed variables dxcC and dxP'^. If we use (6.2) and ( |6.4[ ), then the following inequality 
is easily obtained: 

\dxD^{t,x)-dxD’^{t,x')\ + \dxB’^{t,x)-dxB^{t,x')\ 

< Ct{\x - x'\ + \dxa"-{t,x) - dxa^{t,x')\ + \dxl3'^{t,x) - dxli"'{t,x')\). 


(6.24) 
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Hence, the second result of Lemma 6.8 implies the following: 

|V,X"(s;t,z) - + |V,y”(s;i,z) - V,F"(s; z')l 

<Ct\z-z'\+Ct [ ma^){a,X^{a;t,z))-{d,a”)ia,X^{a-t,z'))\da 

J S 

+ Ct f |(9./3")(a,X"(a;i,z)) - (5,/3")(a, X"(a; i, z'))l 

J S 

Lemma 6.9. Consider the characteristic equations ( |6.21 1 and 

VM = -cosr(r,r(T)), C{t)=x, 


dr 

d 

dr 


ry" (r) = cos (r, ry" (r)), 7y” (t) = a;. 


For any x ^ x', we have the following estimates: 

{i) |^”(s; t, x) - C{s-, t, x')\ + |ry"(s; t, x) - 7y"(s; t, a:')| <Ct\x-x'\. 

(n) \dxCis-,t,x) - d,:,Cis-,t,x')\ 

<Ct\x-x'\+Ct I \dxP'^{T,C{T;t,x)) - da;/3'^{T,C{T;t,x'))\dT. 

J S 

(in) \d,xV''{s]t,x) - da;r]'^{s;t,x')\ 

<Ct\x-x'\+Ct f |aa;a”(r,?y”(r;t,a;)) - i9a;a”(T,ry”(r;i,x'))| rfT, 

J S 

where the constant Ct does not depend on n. 

Proof. The argument of the proof is straightforward and almost the same with that of Lemma 6.8, 
hence we skip the proof. □ 


Lemma 6.10. Consider the iteration functions (6.1)-(6.4|. For any z ^ z', we have 

\x,r+\t,z)-v,r+\t,z’)\ 

< Ct|(V,/“)(Z"(0;Lz)) - (V,/“)(Z"(0;Lz'))| 

+ Ct\z-z'\+Ct f \{dxa^){s,X^{s-,t,z)) - {dxa^){s,X^{s;t,z'))\ds 
Jo 

+ct [ mn{s,x^is-,t,z))-{dxnis,x^{s-,t,z'))\ds. 

Jo 

Proof. We only consider x derivative of since the calculation for v derivatives is almost same. 
Since r+\t,z) = /“(Z"(0; t, z)), we have 

\d,r+\t,z)-d,r+\t,z')\ 

< |(V,/“)(Z"(0;Lz)) - (V,/“)(Z"(0;t,z'))||5,Z”(0;Lz)| 

+ |(V.r”)(Z"(0; t, z'))l t, z) - d^Z^iO; t, z')| 

< Ct|(V./“)(Z"(0;Lz)) - (V./“)(Z"(0;Lz'))| 

+ Ct\z- z'\+Ct [ \{d,xa"-){s,X'^{s;t,z)) - {dxa"){,s,X'^{s;t,z'))\ds 
Jo 

+ Ct f |(a,/3")(s,W"(s;t,z))-(a,/?")(s,X"(s;t,z'))|ds, 

Jo 

where we used Lemma 6.5 (ii), boundedness of norm of initial data, and Remark 6.1. This 
completes the proof. □ 
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Lemma 6.11. Consider the iteration functions (6.1)-(6.4). For any x x', we have the following 
estimates: 

{i) |i9a;a”(t,a;) - 5a;a”(t,x')| 

< CTman{C{0;t,x)) - id,an{C{0;t,x'))\ 

+ Ct\x-x'\+Ct [ \{dj;a^){T,C{T;t,x)) - {dxa^){T,C{T;t,x'))\dT 
Jo 

+ Ct f \{d.nir, r(r; t, x)) - {d^nir, Cir; t, x'))\dT 
Jo 

+ C't/ J \id^DiT,C{T]t,x),v) - id^P){T,C{T]t,x'),v)\dvdT. 

(a) \dj;/3^{t,x) - dx^'^{t,x')\ 

< C'T|(9./3“)(r;"(0;i,x)) - (5,^™)(77"(0;t, x'))| 

+ Ct\x-x'\+Ct [ |(9a;a")(T,?7"(T;t,x)) - (i9a;a”)(T,?7"(T;i,x'))| dr 

+ Ct [ \id^l3‘^)iT,r]‘^{T]t,x)) - {d^j3‘^)iT,r]‘^{T]t,x'))\dT 
Jo 

+ C't/ J \{dj;f^){T,7j^{T;t,x),v) - {dj;f^){T,7j'^{T;t,x'),v)\dvdT. 

Proof. We first consider dxjf‘{t,x) — dxjf{t,x'), i = 0,1,2, where Jq = p", which can be easily 
estimated by (|6.2|) as follows: 


\dxji{t,x) - dxffit,x')\ < Clx- x'l + J\dxPit,x,v) - dxPit^x',v)\dv 

for any i = 0,1, 2. Then, we use (5.14)-(5.16) and the fact that the quantities in (6.11) are uniformly 
bounded by Ct to obtain 

\dxki{t,x) - dxki{t,x')\ < Ct\x - x'\ + J\dxPit,x,v) - dxP{t,x',v)\dv 

for any i = 0,1,2. The proof is now straightforward. We use (5.6|, take x derivative on it, and 
estimate it for x and x' by using Lemma 6.9 and the above inequality together with the fact that 
the quantities in (6.11) are uniformly bounded by Ct, and then we obtain (i). The second estimate 
is obtained by the same argument, and this completes the proof. □ 


Lemma 6 . 12 . Consider the iteration functions (6.1|-(6.4). The sets of derivatives of the iteration 
functions {V^/"}, {dxCx^}, and {dxP^} are equicontinuous. 

Proof. We use Lemma 6.8-6.11. We first define the following quantities: 

e„(t, 5) := sup{|V 2 /”(t, z) - ^zTiJ, z')\:\z- z'\< i5}, 

0n{t,6) := sup{|5xQ:”(t,x) - dxa'^{t,x')\ + \dxl3"^it,x) - dxl3^it,x')\ : \x - x'\ < ^}. 

Note that for any positive integer N we have 

£nit,S) < e„(t, iV(5) < Nen{t,6). 

On the other hand, we can see that for any n 

£„(0,(5) = eo(0,(5) =: eo(<5) and 6>„(0, J) = 6»o(0, ^) =: 0o((5), 
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which are determined by given initial data. By Lemma 6.8 and 6.9, we can choose a positive 
integer M such that 

\Z'^{s]t,z) - Z"(s;t,z')l + - C‘{s\t,x')\ + \rf'{s\t,x) - iif{s]t,x')\ < M6 

for any 6 > 0, \z — z'\ < S, n > 0, and 0 < s < t < T. We now fix a 5 > 0 and use and 

6n{t,6) to rewrite Lemma 6.10 and 6.11 as follows: 

en+iit,S) < CTSoiMS) + CtS + Ct f 9n{s,M5)ds 

Jo 

^ “f Cx^ [ (5) ds 

Jo 

< Cx + £o(i^) + J 9n{s, 6) ds^ , 

and similarly 

0n{t,6) <2CxOo{MS) + 2CxS + 2Cx [ 6»„(r, MJ) dr + 2C't / [ en{T,MS)dvdT 

Jo Jo J\v\<P(t) 

< Ct + do(<5) + J 9n{T,6)dT + j e„(r, d)dr^. 

We apply Gronwall’s inequality to the last inequality and then iterate the above two inequalities to 
conclude that for large n we have 

£n{t, S) + 9n{t, S) < Cx{d + Sq{5) + 9q(S)) 

for any 0 < t < T, and this implies that {V^/"}, {d^a^}, and are equicontinuous. □ 


Lemma 6.13. Consider the iteration functions (6.1 )-( |6.4[ ). Then we have 

iiv,r+i(t)-v,/™+i(t)iioo 

< Enm + Cx f \\dxa^{s) - a 2 ,a'"(s)||oo + Px/^”(s) - i9a;/3"'(s)||oo ds, 

Jo 

where Snm is a small positive quantity which tends to zero as n,m ^ oo and depends on T but not 
on X and v. 

Proof Since /”+^(t,z) = /“(Z"(0; t, z)), we have 

< |(V,/“)(Z-(0;t,z)) - (V,/“)(Z-(0;t,z))||V,Z"(0;t,z)| 

+ |(V,/“)(Z™(0; t, z))| I V,Z"(0; t, z) - V,Z“(0; t,z)\. 

We use Lemma 6.5 (ii), {in), and then ( 6.24[ ) and Lemma 6.7 to obtain 
\VJ^+\t,z)-VJ^+\t,z)\ 

< + C'T|(V,/“)(Z"(0;t,z)) - (V,/“)(Z™(0;t, z))| 

+ Cx [ |(a,a")(r,X"(r))-(9,a")(r,X’-(T))|dr 
+ Cx [ |(a,/3")(T,X"(r))-(9,^")(r,X-(r))|dT 

+ Ct / ||a,a"(T)-a,a"^(T)||oo + ||5,/3"(T)-a.^"^(T)||^ 

^0 
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where we also used Lemma 6.5 (i). We now use the equicontinuity of {V^/”}, {d^a’^}, and 
in Lemma 6.12. Since — Z™ is £nm, the second, third, and fourth quantities in the RHS of the 
above inequality are also Snm which is independent of dxCt’^ and dxP^■ This completes the proof of 
the lemma. □ 


Lemma 6.14. Consider the iteration functions (6.1 )-( |6.4[ ). Then we have 
\\dxa^{t) - dxa^{t)\U + l|5./3"(t) - 9./3'"(t)||oo 

<£nm + CT [ (s) - SajO™(s) 11oo + 11(s) - 5,^/3™ («) I loo ds 

+ Ct [ ||V,r(s)-V,/”^(s)|Uds, 

Jo 

where £nm is a small positive quantity which tends to zero as n,m ^ oo and depends on T but not 
on X and v. 

Proof. As in the proof of Lemma 6.13, we use the equicontinuity of {dxCe'^} and {9x/3"}- We take x 
derivative on (5.6), and then apply Lemma 6.6 to have 

\dxa'^{t,x) - dxa""{t,x)\ 

<£nTn + CT f \\dxoC[s) - dxa^{s)\\ao + \\dxP'^(.s) - dxl3'^(s)\\oods 

Jo 

2 .t 

+ CtY^ / \{dxkf)iT,C{r))-{dxkf){T,r{r))\dr 

i = odO 

+ Ct E f\\(9.kf){T) - {dxkDMlU dr. 

.■_n JO 


i=0 ' 


Since {Vz/"} is equicontinuous, so are {dxkf}, i = 0, 1,2, due to (5.14)-(5.16|. Hence, the second 
integral in the RHS of the above inequality is £nm by Lemma 6.6 (i). Moreover, (5.14)-(5.16) implies 
that the last quantity above is bounded by Snm and ||Vz/”(r) — V2 /'"(t)||oo, and this gives the 
desired result. □ 


Convergence of the derivatives. We first define the following quantities for simplicity: 

frit) :=iiVzr(t)-Vz/™(t)iioo, 

erit) := Wdxa^it) - dxa^{t)r + P./3"(t) - a,/3”^(t)|U- 

Then, Lemma 6.13 and 6.14 are rewritten as follows: 

< £nm + Ct f eris) ds, 

Jo 

ent) < Snm + Ct [ eris) + fris) ds, 

Jo 

which give the following two integral inequalities: 

(6.25) ^(n+l)(m+l)(^) < ^ r 

Jo 

(6.26) erit) < £nm +Ct [ fris) ds. 

JO 
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We fix the constant Ct in (6.25)-(6.26) and iterate (6.25) as in the previous results to obtain 

< £nm + (-'Tt£(n-l)(m-l) H-^-£(n-2)(m-2) + ' ’ ’ H-- £(n-k){m-k) 




^ e ^ maxl^T^m: ■ * * 1 £^(n —fc)(m—fc)} “t” 


(C'Tt)''+^ 


II „(n-fc)(m-fe)|| 

ll/l I|L“[0,T]- 


(fc + 1)! 

Since /{"■ *^^*'"* is bounded by a constant, which depends only on T, and the quantity 
converges to zero as fc —?► oo, the above estimate implies that {V^/”} is Cauchy on [0,T]. Finally, 
(6.26) implies that {dj;a^} and {d^fd^} are also Cauchy, and therefore we conclude that the solution 
we constructed in Section 6.3.1 is C^. 
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